Magnetic field-induced transition in BaVS3 
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The metal-insulator transition (MIT) of BaVS3 is suppressed under pressure and above the crit- 
ical pressure of p cr ~ 2GPa the metallic phase is stabilized. We present the results of detailed 
magnetoresistivity measurements carried out at pressures near the critical value, in magnetic fields 
up to B = 12T. We found that slightly below the critical pressure the structural tetramerization - 
which drives the MIT - is combined with the onset of magnetic correlations. If the zero- field tran- 
sition temperature is suppressed to a sufficiently low value (Tmi <15K), the system can be driven 
into the metallic state by application of magnetic field. The main effect is not the reduction of Tmi 
with increasing B, but rather the broadening of the transition due to the applied magnetic field. 
We tentatively ascribe this phenomenon to the influence on the magnetic structure coupled to the 
bond-order of the tetramers. 

PACS numbers: 71.30.+h,71.27.+a,72.80.Ga 



The nature of the metal-insulator transition (MIT) of 
BaVS3 has presented riddles ever since its discovery a 
quarter of a century agoi^ The ambient pressure MIT 
is spectacular in the susceptibility cusp^* 3 -^ but it is 
also manifested in specific heat, 5 thermal expansion, 3 
resistivity^* 2 -* 3 -^ and optical reflectivity 6,7 anomalies. The 
measurement of transport properties under pressure 2 -* 3 -*^ 
revealed that the insulating phase is bounded by a line 
of MITs in the p-T plane. The transition temperature 
Tmi drops to zero at the critical pressure p cr « 2GPa and 
the high-pressure p > p CY phase is metallic down to the 
lowest accessible temperatures 4 

The order of the MI phase transition was debated, and 
the options of weakly first order, second order, and su- 
percritical transition have all been put forward. How- 
ever, resistivity 8 and in particular magnetoresistivity 9 
measurements under pressure led us to favor a second 
order transition. The character of the transition remains 
unchanged under pressures up to 1.5 GPa as it always 
corresponds to a sharp spike in the logarithmic derivative 
of the resistivity, dlogp/d(l/T). This suggests a line of 
second order MITs. One purpose of the present study to 
find out whether the same phase boundary is preserved 
above 1.5 GPa extending up to the critical pressure. 

A critical line is possible only if it separates phases of 
different symmetry, so the MIT in BaVSs must be at the 
same time an ordering transition breaking some global 
symmetry J£ However, the ambient-pressure MIT at 69K 
is not accompanied by the appearance of detectable mag- 
netic moments, and the first two decades of search failed 
to reveal a symmetry- changing distortion at Tmi-— Fi- 
nally, in 2002 Inami et al. reported the doubling of the 
unit cell along the trigonal c-axis^ 2 - Since the unit cell of 



the T > Tmi structure contains two V atoms along the 
c- axis, 13 the T < Tmi unit cell must contain four, hence 
the symmetry breaking change corresponds to tetramer- 
ization. More recent X-ray scattering measurements re- 
confirmed the existence of four inequivalent V sites at low 
temperatures i 14 i 15 i 16 The X-ray measurements were ex- 
tended to 5K, well below the third and lowest transition 
temperature Tx = 30K of BaVS3. The mysterious Tx- 
transition leads to a low-T phase with incommensurate 
magnetic long range order ji£ and possibly orbital order ji£ 
but the symmetry of the lattice is held at Tx i 14]15]16 

The pressure range 1.5 GPa < p < p cr , where Tmi(p) 
drops below « 20K, had not been studied in detail. It 
is expected that novel behavior may arise if the struc- 
tural transition and the magnetic transition eventually 
meet, or even only get near. As Tx was found to be 
pressure independent in the pressure range investigated 
(p < 1 GPa),— the Tmi(p) and Tx(p) lines may merge, 
which could change the nature of the MIT. 

Another consideration is that even the well-studied 
low-pressure MIT has a dual nature: it is the boundary of 
a phase with spontaneous symmetry breaking (tetramer- 
ization), and at the same time a metal-insulator phase 
boundary. Since the MIT coincides with a symmetry 
breaking transition, the phase boundary may in princi- 
ple be continued to a zero-temperature quantum criti- 
cal point. On the other hand general arguments suggest 
that a metal-insulator transition is likely to become first- 
order before the predicted gap could become too small.— 
Furthermore, it has been pointed out that the quantum 
critical behavior of an MIT may be quite unlike the quan- 
tum criticality of a mere symmetry breaking transition. 21 
BaVS3 has both aspects, thus it presents a situation of 
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FIG. 1: (Color online) Pressure dependence of the metal-insulator 
phase boundary, as identified by the peak in the logarithmic deriva- 
tive of the resistivity. Results for two samples, marked by different 
symbols, are presented. The boxed range indicates the interval 
where results presented in the following figures were obtained. 



general interest but little firm knowledge. 

In this paper we investigate if there is a change in 
the nature of the MIT at pressures close to the criti- 
cal value. In order to get detailed information about the 
range where Tmi(p) is suppressed to zero, we carried out 
magnetoresistance measurements at pressures where the 
characteristic energy of the applied field is comparable 
to that of the transition temperature. Results related to 
the MIT are presented here, while studies about the high- 
pressure metallic phase will be reported in a companion 
papers 

Single crystals of BaVSs were grown by Tellurium flux 
method. The crystals, obtained from the flux by subli- 
mation, have typical dimensions of 3 x 0.5 x 0.5 mm 3 . 
The sample was inserted into a self-clamping cell with 
kerosene as a pressure medium. The pressure was moni- 
tored in-situ by an InSb sensor. The pressure was stable 
within 0.05 kbar in the pressure and temperature ranges 
investigated in this work. The resistivity of the single 
crystal was measured in a standard four probe arrange- 
ment. The current was kept low enough to avoid the 
self-heating of the sample. Magnetic fields up to 12 T 
were applied perpendicular to the current, flowing along 
the crystallographic c direction. 

The zero-field pressure dependence of the metal- 
insulator phase boundary is shown in Fig.[TJ The bound- 
ary consists of two distinguishable parts: Tmi decreases 
with p almost perfectly linearly to about p « 1.8GPa 
where the insulator becomes much softer, and the bound- 
ary begins to drop steeply towards zero, signifying a new 
regime of the MIT. The almost vertical drop of the Tmi(p) 
curve makes it difficult to map this region in detail. The 
points in the boxed range of Fig. 1 represent the data ob- 
tained by attempting to fine-tune the pressure up to p cr . 
Even at p = 1.98GPa where Tmi is as low as to 7.5 K, the 
usual characterization of the transition by d\ogp/d(l/T) 
or magnetoresistivity curves yielded no striking difference 
from the behavior seen up to 1.5GPa. 9 By these crite- 



ria, the transition at 1.98GPa is still second-order. How- 
ever, the considerable smearing out of the phase tran- 
sition by magnetic field - in addition to the downturn 
of the phase boundary - indicates that the character of 
the MIT changes in the parameter range p > 1.8 GPa, 
T < 20K (this domain is not sharply delineated). 

Former magnetotransport studies up to 1.5 GPa 
showed that the low-field phase boundary follows a 
quadratic shaped i.e. ATmi(p, B)/Tmi(p) oc [B / 'B c (p)] 2 , 
where pressure dependent critical field and the zero- 
field transition temperature scales together according to 
HbB c (p) = 1.7- ksTMiip)- In the case of the p = 1.5 GPa 
measurement - and for any lower pressures - no broad- 
ening of the transition was observed up to the high- 
est field investigated B max = 12T. Taking Tmi(p) as 
the characteristic energy scale, this field corresponds to 
5 max « 1/5Tmi- In contrast, in the high-pressure range 
p > 1.8 GPa the transition is considerably smeared out by 
magnetic fields smaller than 1/5 Tmi, as shown in Figs. 2. 
This implies that in the new regime of the transition, i.e. 
in the boxed range of Fig. [TJ the insulating phase became 
considerable less stable against the change of pressure 
and the application of magnetic field. We guess that the 
arising of an electronically soft state of matter is due to 
the presence of competing orders. 

We also found that if Tmi is suppressed below ~ 15K 
then the metallic state can be reached by applying mag- 
netic field. Such field-driven metal-insulator transi- 
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FIG. 2: (Color online) Temperature dependence of the resistance 
measured in presence of various magnetic fields at p = 1.98GPa. 
At this pressure with increasing magnetic fields BaVS3 undergoes 
a field-driven insulator-to-metal transition. 
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tions are shown in Figs. 2 and 3 for p = 1.98GPa and 
p = 1.8GPa, where at zero magnetic field the transition 
temperatures are Tmi = 7.5 K and 12 K, respectively. 
Figure El shows the temperature dependencies at vari- 
ous magnetic fields (T-sweeps at p = 1.98GPa). At low 
fields the low-temperature resistivity is monotonically in- 
creasing with decreasing temperature, indicating an in- 
sulating state. However, at high fields the temperature 
dependence changes character and the ground state of 
the system is metallic. Figure [3] present an example for 
the magnetic field dependence of the resistance at var- 
ious temperatures (B-sweeps at p = 1.8GPa). In both 
cases the field-induced insulator-to-metal transition have 
taken place at values well below B c (p) derived from the 
low-pressure scaling rule. These would be B c = 30T and 
18T for p = 1.8GPa and p = 1.98GPa, respectively. This 
also implies that the transition in the boxed regime of 
Fig. [1] has a different character, and we investigate this 
point in more detail below. 

The overall character of the field-induced transition is 
revealed by the 3-dimensional plot of log R over the B- 
T plane. This is shown for p = 1.8GPa in the upper 
panel of Fig. [H Here the R(B, T) surface is constructed 
from the experimental 5-sweep and T-sweeps; beside the 
results shown in Fig.[3]this plot includes also temperature 
dependencies at various magnetic fields. It is clear that 
at B = the system is most insulating at low T, while 
the drive towards the conducting state looks quite similar 
in B and T directions. 

The transition temperature is generally identified with 
the position of the peak in the logarithmic derivative of 
the resistivity. The middle and bottom panels of Fig. [H 
show the d\og R/d(l/T) curves over the B-T plane. 
The top view of the color-coded surface illustrates the 
magnetic field dependence of the transition temperature. 
Though this might imply that a certain point Tmi (B) 
drops sharply as a function of the magnetic field and a 
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FIG. 3: (Color online) Magnetic field dependence of the resistance 
measured at various temperatures at p = 1.8GPa. (The data are 
normalized to the zero field values.) The low temperature curves 
examples the field-driven insulator-to-metal transition. For com- 
parison the magnetoresistance measured at T = 25K is also shown. 
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FIG. 4: (Color online) 3D plots obtained by temperature and 
magnetic field sweeps at p = 1.8GPa (from light to dark: from 
highest to lowest values). Upper panel: log p over the B-T plane. 
The field gradually suppresses the insulating character. Middle 
panel: Plotting dlogp/d(l/T) over the B-T plane shows that the 
field-induced transition gets weaker and less easily defined at higher 
fields. Middle pan el: Top view of the middle panel, i.e. the magnetic 
field dependence of the transition temperature derived from cusp- 
like anomaly in the temperature dependence of dlog p/d(l/T) 



field-induced first order transition takes placed, the side 
view of the surface reveals that that the main effect of 
the magnetic field is not the shifting of the cusp in the 
derivative, but rather its general suppression. 

The same analysis of the data recorded at p = 1.98GPa 
(Fig. Ej) shows that with increasing pressure the insulat- 
ing domain shrinks, the transition temperature is shifted 
down, but the general features remain unchanged. Again, 
by increasing £?, the transition becomes ill-defined indi- 
cating that the character of the high-pressure MIT is 
definitely different from that of the low-pressure (say 
p < 1.5GPa) sharp transition. We interpret the observed 
broadening as a sign that magnetic ordering phenomena 
(related to those appearing at the low-pressure ordering 
transition at Tx) begin to interfere with the MIT, from 
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FIG. 5: (Color online) 3D plots obtained by temperature and 
magnetic field sweeps at p = 1.98GPa. Upper panel: log p over the 
B - T plane, Lower panels: d\og p/d(l/T) over the B - T plane 
(from light to dark: from highest to lowest values). A comparison 
to Fig. 4 reveals that with increasing pressure the insulating domain 
shrinks. 



p ~ 1.8GPa upwards or Tmi ~ 20K downwards. 

The suppression of non-ferromagnetic order by mag- 
netic field is a common phenomenon but it is not neces- 
sarily associated with the broadening of the transition. 
There are many systems (including ordinary antiferro- 
magnets) which possess breakable symmetries even in the 
presence of an external magnetic field. This residual sym- 
metry can be broken by an order parameter which is not 
induced by the field. The critical temperature decreases 
with increasing B but the transition remains well-defined. 
An example is given by the spin density wave system 
of (TMTSF)2PFq which continues to undergo a sharp 
phase transition around T =10K even in magnetic fields 
as high as B w 30T^ 

BaVS3 is definitely not in the above category, the 
metal-insulator transition is considerably broadened in 
fields as demonstrated in Figs. 4 and 5. There are two 
possibilities: either the order of the phase is induced by 



the field (thus in the presence of the external field there 
is no symmetry breaking and the transition becomes ill- 
defined), or the observed phenomena are kinetic, and can- 
not be explained in the language of equilibrium phase 
transitions. 

Taking the first possibility, field induces in linear or- 
der the corresponding components of spin and orbital 
momentum. This is relevant if the low-T phase BaVS3 
possesses a net ferromagnetic moment (since the ferro- 
magnetic moment may be very weak, it may have es- 
caped experimental detection). Even if the zero-field 
ground state carries no total magnetization, field could 
also induce higher-order uniform moments. If one these 
is mixed into the order parameter of BaVS3 at T < Tx, 
the observed broadening may once again follow. 

Future experiments will decide whether the broadening 
of the transition can be explained in terms of equilibrium 
phases, as described above. However, we wish to point it 
out that there is an alternative: namely that the broad- 
ening is a relaxation effect, and thus not explicable within 
the framework of equilibrium phase transitions in ideal 
systems. 

There is an independent observation showing that the 
character of the MIT changes, and becomes mixed with, 
or preceeded by, magnetic ordering phenomena in this p- 
T regime; above p ~ 1.8GPa hysteresis loops may appear 
both in £>-cycles and in T-cycles2*22. The shape of the 
hysteresis loops depends on the sweep rate, suggesting 
that we are dealing with a relaxation phenomenon. Hys- 
teresis appears both above and below the critical pres- 
sure, thus it is plausible to relate it to the complicated 
long-period order which was observed by magnetic neu- 
tron scatterin g 17 ! 18 ! 19 . Similar, so called weak relaxation 
effects were observed in incommensurate density wave 
systems and successfully interpreted in the picture of a 
non-ideal system going through a hierarchy of pinned 
states^. Formally, the random impurity pinning leads 
to a broad distribution of the relaxation times, and this 
results in logarithmic, power law, or stretched exponen- 
tial relaxation. The finding that the ambient-pressure 
low-T phase of BaVSs displays a magnetic order which 
is incommensurate both in the basal plane and along the 
c-axis^ may provide the key elements for an analogous 
relaxation mechanism. If an incommensurate structure is 
locally pinned to defects, while the optimal wave number 
is temperature dependent, the temperature sweep yields 
an inhomogeneous state which try to relax towards the 
optimal structure with a locally varying rate. 

To summarize, we investigated the high-pressure low- 
temperature region where the insulating phases are sup- 
pressed. The zero field boundary was followed down to 
Tmi ~ 8K. We did not directly see a first order jump 
in either the temperature, or the pressure, or the field 
dependence, but in the range Tmi < 20K, p > 1.8GPa 
we observed a considerable smearing out of the phase 
transition by magnetic field. We suggest this is a sign 
that the metal-insulator transition and magnetic order- 
ing phenomena appear combined in this (p, T) regime. 
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